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^ '■ Abstract. 

\ Results of 2D numerical simulation of magnetorotational mechanism 

^ ■ of supernova explosion are presented. It is shown that due to the differen- 

•y^ I tial rotation of the star toroidal component of magnetic field appears and 

■ grows with time. Angular momentum transfers outwards by the toroidal 
i magnetic field. With the evolution of the process part of the envelope of 

■ the star is throwing away. The amount of thrown away mass and energy 
, are estimated. The results of the simulation are qualitatively correspond 
' to supernova explosion picture. 

o ■ 

|H '. !• Introduction 

' One of the possible mechanisms of supernova explosions is magnetorotational 

J> . mechanism, suggested by Bisnovatyi-Kogan (1970). It is based on the transfor- 

mation of the gravitational energy to the energy of explosion by magnetic field. 
Poloidal magnetic field in differentially rotating star is twisted and toroidal com- 
I ponent of the magnetic field appears and grows with the time. When the force 

produced by the magnetic field becomes comparable to the gravitational force 
it pushes the matter of the star outward. 

A numerical simulation of the magnetorotaional supernova explosion in ID 
approximation has been made by Bisnovatyi-Kogan et al. (1976), Ardeljan et 
al. (1979), Miiller and Hillebrandt (1979). 

Pioneer 2D simulations of supernova explosion with magnetic field has been 
made by Le Blank and Wilson (1970). They used unrealistically large magnetic 
fields: initial magnetic energy was comparable to the gravitational energy of 
the star. More realistic case was considered by Ohnishi (1985). As result of the 
simulations they got axial jets. Simulation of the collapse of magnetized rotating 
gas cloud in 2D has been done in number of papers, see for reference Ardeljan 
et al. (1996b). 

Realistic magnetic fields in the stars are rather small (the ratio between 
magnetic and gravitational energies is about 10~^ — 10"*^. Such a small magnetic 
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fields are the main difficulty for the numerical simulations, because of very small 
hydrodynamic time scale and huge time scale of the magnetic field amplification. 

Numerical simulation of such stiff problem requires application of the im- 
plicit numerical methods, which are free from Courant restriction on the time 
step. 



2. Mathematical formulation of the problem, numerical method 

The problem of the collapse of rotating magnetized gas cloud is described by 
the set of magnetohydrodynamical equations with selfgravitation for the gas 
with infinite conductivity (Landay, Lifshitz (1984)). Axial = 0) symmetry 

and symmetry to the equatorial plane are assumed. Here ^ = ^ + u • V , 
X = {r,ip, z) , u - velocity vector, p - density, p - pressure, H - magnetic field 
vector, $ - gravitational potential, e - internal energy, G - gravitational constant, 
TZ - universal gas constant, 7 - adiabatic index, H (X" H — tensor of the range 2. 
The problem was solved in the restricted domain. Outside the domain a density 
was taken close to zero. 

For the initial conditions, we assume that the cloud is rigidly rotating uni- 
form gas sphere with the following parameters: 

r = 3.81 • lO^^cm, p = 1.492 • lO'^'^g/cm^, 
M = 1.73M0 = 3.457 • lO^^g, 7 = 5/3, u"^ = = 0, 

Erot0/\Egr0\ = 0.04, ^ino/l^grol = 0.01, ^^magl / 1 ^inl | = 0.05. 

The nondimensional and scaling values were chosen as: 

r = fxo, z = zxq, u = miQ, 
p = ppo, e = ££o, T = TTo, $ = i$o, 

no = JiTrGpoxl, to = —, Po = PqUq = Poxlto^, To = 

* no /v 

£0 = -"o = a^o^o ^> Ho = Vpo = xoto^pl^^, ^0 = ATrGpox^. 

Here the values with the index zero are the scale factors and the functions under 
sign "tilde" are dimensionless functions. All plots and figures presented below 
are in nondimensional form. 

Boundary values of the gravitational potential $ have been calculated by 
the formula of volume potential (Courantfe Hilbert, 1962). At r = we assume: 

Vr = = Hr = = curl^H = curl^H = 

and at z = we assume: 

acurl^H 

Vz = Hz = — = 0. 

oz 

For our simulations wc used implicit conservative Lagrangian operator dif- 
ference scheme on triangular grid with reconstruction (Ardeljan et al, 1996a). 
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Figure 1. Triangular grid and density distribution at t = 5.153. 



The theory of this numerical technique (investigation of its approximation, sta- 
bility and convergency) has been developed by Ardeljan et al. (1987). 

Grid restructuring procedure allows us to overcome grid overlapping sit- 
uation in gas flows with nonuniform contractions or expansions or flows with 
vortexes. Due to the implicitness of the scheme it was possible to time steps, 
which are much larger then Courant time step. 

The details of the scheme used in this paper will be published elsewhere. 

3. Results 

Initial number of grid knots was 5000. At the initial time moment (t = 0) 
magnetic field is "switched off', because it does not influence significantly on 
the hydrodynamical collapse stage, rather short in comparison with time of the 
amplification of the magnetic field. 

After few consequent contractions and expansions the cloud comes to the 
differentially rotating stationary state, at ti = 5.153 and consists of rapidly 
rotating dense core and extended light envelope, (Figure The angular velocity 
of the cloud at this time presents a function with maximum of gradient at the 
transition zone between core and envelope (Figure ^). 

At the moment ti = 5.153 initial magnetic field of the following configura- 
tion (Ardeljan et al. 1996b) (Figure. |): 

Hro = Fr{0.5r, 0.5z - 2.5) - Fr{0.5r, 0.5z + 2.5), H^o = 0, 
H^o = F^{0.5r,0.5z - 2.5) - F;,(0.5r, 0.5z + 2.5), 

^r[r,z) - (^(^2 + 1)3 - (^2 + 1)5 j '^^^^'^^ " |^(^2 + l)2 " (^2 + 1)4 j ' 

has been "switched on" with coefficient k = 0.43. This magnetic field is sym- 
metrical about the equatorial plane, it has divHo = 0, but it is not force-free. 
The choice of initial magnetic field configuration was connected with an at- 
tempt to avoid a singularity, which could create big numerical problems. After 
the moment of "switching on" the magnetic field, the matter starts to move to 
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Figure 2. Distribution of the angular velocity along r-axis and inital 
magnetic field at t = 5.153. 



Figure 3. Toroidal magnetic field and velocity field at t = 10.234. 



the periphery of the envelope due to the magnetic force produced by the ini- 
tial configuration of the magnetic field, what leads to appearance of a "finger" 
which grows at the outer boundary of the cloud. This artificial "finger" has 
no significant influence on the evolution of the inner parts of the cloud. Due 
to quadrupole type of symmetry of the poloidal magnetic field the appearing 
toroidal component has 2 extremums. First is situated at the equatorial plane 
in the zone between core and envelope, and the second is situated above the 
equatorial plane close to the z-axis (Figure |3|) . 

At t = 10.234 toroidal part of magnetic energy reaches its maximum. While 
the total energy of the toroidal component of the magnetic field of the cloud is 
much smaller than its gravitational energy, the density of the toroidal magnetic 
energy in the regions near extremums of the becomes comparable to the 
density of the internal energy what pushes out the matter of the cloud. Starting 
from t = 11.342 the kinetic energy of the pushed part of the envelope near 
equatorial plane becomes higher than its potential energy and can be ejected 
from the cloud. The ejection along equatorial plane leads to change of the shape 
of the cloud, which is stretched along r-axis (in r, z coordinates), see (Figure 
related to the last computational moment t = 32.634. 
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Figure 4. Velocity field at t = 32.624. 



Figure 5. Time evolution of the ejected mass in % to the total mass 
of the cloud. 

At Figure ^ time evolution of the ejected mass of the envelope in % to the 
total mass of the cloud is presented. The ejected part of the matter contains 
about 0.09% of the gravitational energy of the stationary cloud after collapse. 

The problem of the collapse of the rotating magnetized gas cloud, described 
above can not be considered as complete and final explanation of the mechanism 
of the supernova explosions. However the results of such simplified formulation 
give evidence in favour of the magnetorotational mechanism can as an explana- 
tion for the problem of the supernova explosions. 

Our recent simulations of the collapse of the star with a real equations of 
state and neutrino loses as in the paper by Ardeljan et al, 1987a show that 
the star contracts much stronger and rotates much more differentially, than it 
was for the simple equation of state used in this paper. We may expect a more 
efficient magnetorotational explosion. 
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